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The critical behavior of a quenched random hypercubic sample of linear size 
L is considered, within the "random- Tc" field-theoretical model, by using the 
renormalization group method. A finite-size scaling behavior is established 
and analyzed near the upper critical dimension d = 4 — e and some universal 

C^ . results are obtained. The problem of self-averaging is clarified for different 
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I. INTRODUCTION 

The description of effects of disorder on the critical behavior of finite-size systems have 
attracted a lot of interest [|l]-§]. Up to now a discussion takes place of whether the introduced 
disorder influences the finite-size scaling (FSS) results 0,|], compared to the standard FSS 
results, known for pure systems P-plH- A formulation of general FSS concepts for the 
case of disorder is strongly complicated due to the additional averaging over the different 
random samples. For a random sample with volume L'^, where L is a linear dimension, 
any observable property X, singular in the thermodynamic limit, has different values for 
different realizations of the randomness and can be considered as stochastic variable with 
mean X and variance (AX)^ := X^ — X , where the over line indicates an average over 
all realizations of the randomness. Here, an important theoretical problem of interest is 
related with the property of self-averaging (SA) |jl2|]. If the system does not exhibit SA a 



measurement performed on a single sample does not give a meaningful result and must be 
repeated on many samples. A numerical study of such a system also will be quite difficult. 
This point has been studied recently by means of FSS arguments [|l],§], renormalization 
group (RG) analysis [Q,|| and Monte Carlo simulations MM- The quantity under inspection 
is the relative variance Rx{L) := {AX^/X . A system is said to exhibit "strong SA" if 
Rx{L) ~ L""^ as L ^ cxd. This is the case if the system is away from criticality, i.e. if 
L ^ ^. At the criticality, i.e. when L <C ^, the situation depends whether the randomness 
is irrelevant {impure > 2/(i, e.g. "pure", P-case) or relevant {I'random > 2/(i, e.g. "random" 
R-case) 0. One calls the former case "weak SA", since Rx{L) ~ /,("/'^)p"'<= ^ and the latter 
case "no SA", since Rx{L) is fixed nonzero universal quantity even in the thermodynamic 
limit 0. 

The lack of SA in disordered systems implies that the standard FSS breaks down. Thus, 
one needs to formulate a FSS theory, suitable for the case under consideration. There is an 
ongoing activity in this field but an understanding of the problems at a deeper level is still 
desirable [|l],|,i. 

The successful application of the field theoretical methods U3|,|H| for the analytic calcu- 



lations of the size- dependent universal scaling functions for pure systems makes the possible 
extension of these methods very appealing also in the case of disordered systems. The ap- 
propriate field theory is the N-component ip"^ theory with a "random- Tc" term |[T5|-pT|. 
The analytic difficulties raised by the disorder usually are avoided with the replica trick by 
solving an effective pure problem [0. However, the perturbative structure of the theory 



is still much more complicated than for the corresponding pure system, arising additional 
difficulties in the applicability of the ideas proposed in |p!3| , p^ . Some of them are of rather 
hard computational nature. 

Between the generic models for magnetism, the most popular and relatively well studied 
by means of Monte Carlo simulations is the Ising model (i.e. N = 1) due to the fact 
that in three dimensions it is the model for which in accordance with the Harris criterion, 
randomness is relevant. But its RG analysis is complicated as a result of the well known 
accidental degeneracy in the recursion relations that needs higher order in e and so the 
use of the loop expansion to second order (see e.g. [|T5ip!8H 20[| ) . This leads to the apparent 



computational difficulties in the finite-size treatment of the system. The e- calculations based 



on the loop expansion to second order are not done even for the pure finite-size systems. In 
this situation, more attractive for the finite-size RG study are the disordered XY (A^ = 2) 
and Heisenberg (A^ = 3) cases as simpler and having the same qualitative features governed 
in d = 4 — e dimensions by a random fixed point. 

Other problems have more basic nature and are related with the breaking of the replica 
symmetry (see [^). Since its deeper understanding is still lacking, they are beyond of our 
interest in the present study. 

In this paper we analyze the finite-size properties of an N-component {N > 1) model of 
randomly diluted magnet with hypercubic geometry of linear size L. Exact calculations are 
performed in the mean-field regime d > 4, and up to the first order in e near the upper critical 
dimension dc = 4. Although in this case, the problem of usefulness of the corresponding 
series expansions away from the dimension d = 4 — e arises and is questionable (see ||2^ 



and refs. therein), we shall show that many generic FSS properties of the model can be 
established and we hope this would have implications for more realistic cases. 

The paper is organized as follow. In Sections |I| and |ITl| we define the model and the ef- 
fective Hamiltonian. In Section ^ we perform the analysis in the zero-mode approximation. 
In Sections and we give the expressions for the shift of the critical temperature and the 
renormalized coupling constants in first order in e. Section |V11| deals with the verification of 



the FSS and the analysis of the problem of Self-Averaging is given in Section |VI1I|. Finally 



in Section [X we present our main conclusions. 



II. MODEL 

We consider the "random - T^" Ginzburg-Landau- Wilson model of disordered ferromag- 
nets (see, e.g. |p!5|-pO[1) 



Hr = -\ j^^ d'x[t\i,{x)\' + v.(=r)|V(=r)r + c\Vi^{x)? + ^mx)\% (2.1) 

where il){x) is a N-component field with ilP{x) = Z^ili i'K^) ^^^ ^^^ random variable (/^(cc) 
has a Gaussian distribution 

with mean 



(fix) = (2.3) 

and variance 



<^{x)<^{x/) = Ar (a? - x') . (2.4) 

The over line in ( p.4| ) indicates a random average performed with the distribution P{ip{x)). 
Here we will consider a system in a finite cube of volume L'^ with periodic boundary condi- 
tions. This means that the following expansion takes place 



^(^) = 7^II^(^)exp(ifc.a:;) 



L'^ 



(2.5) 



and 



^{x) 



L^ 



y^ V^(fc) exp(ifc.a;) 



(2.6) 



where fc is a discrete vector with components ki = 2TTni/L , Ui = 0, ±1, ±2, ..., i = 1, ...,d 
and a cutoff A ~ a~^ (a is the lattice spacing). In this paper, we are interested in the 
continuum hmit, i.e. a — ;> 0. 

In our case of quenched randomness one must average the logarithm of the partition 
function over the Gaussian distribution ([2. 2]) to produce the free energy 



J^lHr 



D<^{x)P{<^{x))lnZr 



where 



Zr = Tr^ exp [Hr 



(2.7) 



(2.8) 



It is well known that the direct average of Tir over the Gaussian leads to equivalent 
results [E4| for the critical behavior as the n = limit of the following "pure" translationally 



invariant model 25 



Hp{n) 



1 ^ 



d'x\M^)\'\M^)\']. (2.9) 



A " 

+- H 

Here ipaix), a = 1, . . . ,n [n beings the number of replicas) are components of an {n x A^)- 
components field il'ix). Because of this equivalence, the model Tip has been the object of 
intensive field-theoretical studies (see [^] and refs. therein) in the bulk case. Much less is 
known for the equivalence of Tir and the n = 0-limit of Tip in the finite - size case. Problems 
may arise when finite-size techniques are used, since both procedures L ^ oo and removing 
of disorder by the "trick" ra — > may not commute. 



III. THE EFFECTIVE HAMILTONIAN 



In this work we will use the RG technique introduced in [|T^] and [Q for studying pure 
systems with finite geometry. This technique permits explicit analytical calculations above 
and in the neighborhood of the upper critical dimension. The main idea is to expand the 
field in (|2.1|) in Fourier modes and then to treat the zero mode separately from the nonzero 
modes. The nonzero modes can be treated by the methods developed for the bulk systems 
(e.g. loop expansion), while the zero mode, whose fluctuations are damped at the critical 
temperature, has to be treated exactly. This overcomes the problems due to the infra red 
(IR) divergences that take place in finite size systems (see reference 



In our more complicated case we have two possibilities: to consider the random model 
Eq. (|2.1| ) or to consider the replicated pure model Eq. (^]9|). The last one is closer to 
the case treated in |]13|] and [|I^ by getting around the difficulties due to random average 
performed with P{ip{x)) and is used in the present study. For this case, the replicated 
partition function is given by 



Zp{n) 



V'ipexp ['Hp{n)] . 



(3.1) 



We decompose the field ipi^) i^ito a zero momentum component = L^'^ J d'^xipix), which 
plays the role of the uniform magnetization and a second part depending upon the non-zero 
modes a = L~'^X]fc^o'^(^) 6xp(— ifc.a?). After some algebra, the partition function can be 
expressed as 



Zp{n) 



T>(j)V(7 exp 



L 



d n 



E{r«A + ^ji) + ^(±^i 



a=l 



12^° 



\a=l 



1 " 



ro + fc^ + 



2 , ""0 , 2 ^0 -A ,2 



E' 

/3=1 



cr^ + higher powers of a 



(3.2) 



Here the terms involving / d'^xa vanishes since a depends only on non zero modes. The terms 
containing a are treated using diagram expansion, leading to the effective Hamiltonian in 
the one-loop approximation 



Td n / 



u[n) 



2 ^ \ -■--" ■ 12 ^" 

a=l \ / 



+ 



L'^A(n) 



E' 



^^V, 



(3.3) 



where t{n),u{n) and A(n) will be presented bellow. With the help of the identity 

faA^\ 1 



exp 



(27ra)V2 



dyexp[-{l/2a)y^ + yA], 



we get 



/oo r /"oo 

_JyPn{y) [SnI ci|0| 101^-1 exp(Hf(n)) 



where 



nf{n) = --L'^ 



i{n) 



ld/2 



' 12 ^ 



(3.4) 



(3.5) 



(3.6) 



is an effective Hamiltonian with a random variable y with Gaussian distribution (depending 
on A{n)) 



Pn{y) 



exp 



2A(n) 



\j2'Kl\{n) 



(3.7) 



and S 



N 



r(Ar/2) 



is the surface of a A^-dimensional unit sphere. 



Let us note that the above mentioned equivalence between the models (|2.1| ) and (|2.9| ) in 
the used approximation may be mathematically expressed, within the used approximation 
by the following relation: 



^[^r] 



On ^ 



Znin) \n=0 



From Eqs. ( p. 51) and (|37^ ), and by using the identity 



|-A"(n)|„=o=lnA(0) 
on 



for the free energy, we get 



J^lHr 



dyPoiy)\nZriO), 



where 



/•oo 
Zr{0)=SN rf|0||0r-^Xp[Hf(O)] 

Jo 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



is the partition function for the random system (|3.6|) after taking the limit n ^ 0. The 
obtained effective "random-Tc model" ( p.6|) , distributed with Gaussian weight ( p.7|) , is the 
analytic basis of this paper. The effective constants t{n), u{n) and A(n) involve n and finite- 
size L as parameters. For describing the finite-size properties of the initial model ( p.l| ), as 
follows from Eqs. ( |3.10 ) and ( 3.11 ), it is necessary to set n to zero. In the next sections we 
shall consider the results of this procedure. 



IV. THE FSS EXPRESSION FOR THE FREE ENERGY AND CUMULANTS IN 

THE ZERO MODE APPROXIMATION 



If we neglect the loop corrections this corresponds to the mean-field approximation. Then 
the zero mode playing the role of the uniform magnetization may be treated exactly. In this 
case the effective Hamiltonian (13.31) of the model reduces to 



n, =--L 



n -in A / '^ ^ 



(4.1) 



Now using an appropriate rescaling of the field |0| = (mL'^) ^/'*$ and introducing the 
scaling variable 

/i = tL''l''u~^!\ (4.2) 

for the partition function Eq. ( |3.11|) in the mean-field approximation we obtain 

Zf^{^) = {uL'^r^/'X^ii, + ylu'1% (4.3) 



where we have introduced the following auxiliary function 



In(z) 



POO ( 1 

S^\ rf<l><l>(^-i)exp -- 
JO I 2 

From Eqs. (|3.10|) and ( [4.3|) we get for the free energy 
1 



^$2 _^ $ 

12 



j^ 



U 



MF 



If we introduce a second scaling variable 



c/yexp (-2^) l'^ 



(mL'^)-^/^X^(/x + 2//w^ 



/2 



u 
and using Eqs. (|A4|) (see appendix), Eq. (|4.5| ) takes its final form 

1 



(4.4) 



(4.5) 



(4.6) 



J^ 



n 



MF 



V2^ 



rf3,e-(-/0V(2A)i^ 



D_N„{V^x)\--^{\ + ^?) + ^\n( 



12^ 



(4.7) 



For the Ising case A^ = 1 a similar expression for the quenched free energy in a slightly 



different context is obtained and its analytic structure is studied in references [^^HUl- Ob- 
viously ( |4.7| ) is well defined for any positive A and in the limit A ^ we recover the well 
known result for the free energy of the pure model. 

In addition to the free energy, one also needs to know the correlation functions. Within 



the replica method the averages of the fields {0/?} are defined by (see e.g. |T2 



(10/3^™) 



n} 



where 



lim 



tMF/ 



Zi 



'N 



n 

\a=l 



{\<i>a\r-\\My"'exp{n; 



Zp' (n) 



S-^ J (j[ d\<P^\\ {\<P^\f-' exp{H 



MF\ 
P ' 



(4.8) 



(4.9) 



Note that the final result must be independent of replica index /3, because Ti^^^ is invariant 
under permutation of the replicas. After taking the limit n — ;> 0, we end with the following 
expression: 



M2m := (|0/3p™)«MF 

In a similar way 



V2^ 



Xm(x) 



(4.10) 



(-M2) 



(I0«n0/3P) 



{uL 



d\-l 



Hi. 



V2^ J- 



dx 



J, 



JV+2W 
In{x) 



-(x-/.)V2A 



(4.11) 



From Eqs. ( [1.10| ) and ([4.11|) , when /x = and iV = 1 we obtain the results of Ref. 



In terms of the normalized magnetization Ai the susceptibihty is given as 

X = L''M~2. (4.12) 

Another quantities of importance for numerical analysis is the Binder cumulant defined by 



S^l-i^ (4.13) 

and the cumulant, specific for the random system defined as 

Since the parameter R is the relative variance of the observable (the susceptibility), as 
we said in the Introduction, it is a measure of the self-averaging in the random system. If 
self- averaging takes place this quantity should be zero in the thermodynamic limit. 



V. FINITE-SIZE SHIFT OF T^ : LOWEST ORDER IN e 

The perturbatively calculated parts of the free energy and cumulants, which contain 
contributions of all nonzero modes, depend, to one-loop order, on the shift of the critical 
temperature and on the renormalized coupling constants u and A. The application of the 
finite-size e-expansion to the model system ( |2.9| ) requires the corresponding renormalization 
constants. 

To one loop order, using the minimal subtraction scheme, before taking the n — > limit, 
we obtain: 

Z,.l,^,_i±!!^A, (5.1a, 

Z„^1^4^u-'A (5.1b) 

be e 

N + 2^ 8 + nN ^^ ,_ , 

Za = 1 + ^^u A. (5.1c) 

6e It 



InEqs. dp 



■(47r)^/2r(rf/2)' 
The /3 functions associated to u and A have the form 

iV + 8.2 



^ ~ ''(47r)'^/2r(rf/2)' ^^-^^^ 

A = A— —J———. (5.2b) 



/3„ = —ue H — u — 6-uA, (5.3a) 

6 

;9a = -Ae - —- A^ + ^—ii^. (5.3b) 

8 



The fixed points of tliis system first have been studied in pOi and for the purposes of the 



impurity problem in |17|. The values of m and A in the fixed point, interesting in the random 
case, are 



u [n] 



/\*{n) 



6(4 - nN) 



IQ{N -l)-nN{N + S) 
2{A-N) 



16(iV- l)-nN{N + S) ■ 
The corresponding expression for the exponent v up to the first order of e, is: 

1 6iV(l-n) 



(5.4a) 
(5.4b) 



v{n) 



16(iV 



nN{N + 8) 



e. 



(5.5) 



It should be noted here that we shall consider A^-component fields with 1 < iV < Nc{d), 
where Nc{d) = 4 — 4e + 0{e^) is the critical number of spin components that defines the 
stability of the random fixed point in the n = limit. The stability of the different fixed 



points of the model has been also considered in [^. The analysis of the Ising case {N = 1) 
needs to perform a loop expansion to second order (see Introduction) and is beyond the 
scope of the present study. 

As it was explained above, the loop corrections will be treated perturbatively on the 
nonzero k modes. In the lowest order in e, this procedure generates a shift of the critical 
temperature t -^ t{n); 



i{n)=tZt + tL, 
where the term tZt is coming from the one- loop counterterm (see ( ^.1|) ), and 



N + 2^ 2 + Nn 



6 



-u 



A 



1 
1? 



E' 



1 



fc2+t 



(5.6) 



(5.7) 



is the finite-size correction. The two diagrammatic contributions for t^ are shown in FIG. 
|I|. Both diagrams from the u and the A contributions differ only by their numerical factors. 
The prime in the c?-fold sum in the above equations denotes that the term with a zero 
summation index has been omitted. 

After some algebra (details for the pure case (A = 0) see in flSl), near the upper critical 
dimension d = 4 — e, we obtain 



t{n) = t + 



A^ + 2 . 2 + Nn 



12 



-u — 



A 



(tlnt + 4L-2F4,2(tL=^ 



(5.8) 



where 



Fd,2{x) = / dzexp 
Jo 



xz 



(2vr) 




2^2 



1 - '^''' 



(5.9) 



Some particular values of the constant -Fd,2(0) and a method of calculation are given in |R2 



At the fixed point u = u*{n), A = A*{n), up to the first order in e, the terms proportional 
to In L cancel and Eq. ( ^.8| ) can be written in the following scaling form 



2 .. 3(71 -1)N 



t(n)L' = y \ ^- \y\ny + 4^4 2(2/)! e, (5.10) 

where the scaling variable y = tL^^'^^^'^ has been introduced. 

In the n = limit, the expression for the exponent measuring the divergence of the 



correlation length is |16,17 



UR u{0) 8{N-1 

instead of the critical exponent for the pure case 



^ ^ - 3^ ,5.11) 



- = 2-^. (5,12) 



The scaling form ( ^.10|) can be written for this case as 



S7V 
mL' = y + ^g^^_^^ [ylny + 4F,M] e, (5.13) 

where y = tL^/'^^. 

From the above expression one can obtain the large-L asymptotic form of the Tc(cxo) 
shift, i.e. 

T,(L)-T,(oo)~L-i/'^«. (5.14) 



In Eq. ( |5.14D , Tc(L) := Tc{(p,L) denotes the average pseudo critical temperature {Tc{(p, L) 



is pseudo critical temperature for a specific random realization ^{x)) and Tc(oo) 



lim2,^oo^c(-^)- Eq- (|5.14|) was suggested in p[. Combined with the fenomenological FSS 
theory it gave rise to the lack of SA, and is confirmed by numerical studies (see [^). Here 
it is verified independently and directly. 

VI. RENORMALIZATION OF THE COUPLING CONSTANTS : LOWEST 

ORDER IN e 

We perform the renormalization in a similar way also for «(«) and A{n) by taking into 
account the diagrammatic contributions, from u and A shown in FIG. 0. The result is 

u{n)=uZu + UL (6.1a) 

A{n) = AZA + AL, (6.1b) 

where uZu and AZ^ are the one-loop counterterms for the coupling constants, and 



10 



Ul 



A, 



' 2 ^ + 8 . . 
u — dmA 



uA 



6 



A^ + 2 8 + nN 



1 y^/ 1 



(P + 1) 



2' 



(6.2a) 
(6.2b) 



are the corresponding finite-size corrections. As one can see the summand in Eq. ( |6.2| ) can 
be expressed as the first derivative of the summand of Eq. (|5.7|) with respect to t. So, at 
the fixed point, algebraic transformations similar to those performed in the previous section 
lead to 

hl + \ny)e + 2eF'{y) 



u\n)U = u*{n) 



A*{n]U = A*(n) 



1 



:i + \ny)e + 2eF'{y) 



(6.3a) 
(6.3b) 



where the prime indicates that we have derivative of the function -^4,2(2/) with respect to its 
argument. 

The results for the disordered system simply follow by setting n = 0. From the results for 
the shift of the critical temperature ( p.lOj ) and the renormalization of the coupling constant 
u, given by Eq. ( |6.3aD , we reproduce the results for the pure FSS case, by setting A = 
and n = 0. Moreover, this result still holds even if we find the FSS corrections after the 
analytical continuation to n = 0, expressing the commutativity of the problem. 



VII. VERIFICATION OF FSS 



Let us consider the scaling variables 



^(n) = i{n)L^/y^/ii^, X{n) = A{n)/u{n). 

At the fixed point they can be expressed in terms of scaling variable y = tL^/'^'^") : 



fx^in) 



1 f 1 



u*{n) 



12N{n- 1) 



and 



16{N -l)-nN{N + 8) 

A-N 



{4:-N){4:-nN) 
lQ{N-l)-nN{N + 8) 

F,M^-yF'iy)e 



\ny 



X^in) 



3{A-nN)' 



(7.1) 



(7.2) 



(7.3) 



In the limit n = 0, Eqs. ( |7.2| ) and ( |7.3| ) yield the following scaling variables describing the 
disordered system (^?T|) : 



/x'^ := fi%0) 



1 f 1 

y^i^y 



+ 



u*{0) I 

3N 



4-N , 

1 H ; m y 

4(A^-1) _ 



-F,,2iy)e-yF'{y)e 



(7.4) 



11 



where y = tL^/'^^, and 



4- A^ 
A* := A^(0) = -^. (7.5) 



These equations verify the finite-size scahng hypotheses and show that we are really dealing 
with a one-variable problem, since the second variable A* is a fixed universal number. At 



the critical point t = 0, since the constant -^4,2(0) = —8 In 2, see Ref. |^, we have 



A^ln2 3e ,„ , 

Numerical values for different thermodynamic quantities can be obtained with the help of 
Eqs. ( [7. 51 ) and (|7.(j|). Note that the scaling variable fiQ is proportional to ^/e. Consequently 
all the e-expansion results will be expressed in power of ^/e as it was the case for the pure 



systems (see [jT3[ for example). 



VIII. CUMULANTS AND SELF- AVERAGING 



In Ref. 01, the Binder cumulant B and the the relative variance R (Eqs. ( |4.13| ) and ( |4.14| )) 



are calculated analytically and numerically at the critical point T = Tc in the asymptotic 



regime A ~ 1/4 for A^ = 1 and d = A. In ref. ||3^ (see also |^), the same quantities were 
calculated numerically for A^ = 1 and d = 3. In both cases the results showing that the 
system exhibits a lack of SA. 

In the remainder of this section we concentrate on the calculation of the cumulants B 
and R ( |4.13| ) and ( |4.14| ) in cases d > 4 and d = A — e. The meaning to consider the case d > A 
is in its simple analytical nonperturbative treatment. Although the results based on the e - 
expansion give only a qualitative description of the three dimensional physics, we hope that 
they shed some light at least on the applicability of the theory for studying diluted models. 

Let us first note that if 1 < A^ < 4 and d = 4 — e, the case under consideration applies to 
the situation (see Eq. ( |5.11|) ) where ur > 2/d and randomness is relevant (R-case). Up to 
the first order in e, due to the RG arguments, no SA must be expected near the critical point 
0. This statement is supported also by our RG calculations. In Table | and Table || we 
present the corresponding universal numbers for B and R at rf > 4 and d = 3 in the region 
Lf^^ = J <^ 1, i.e. in the vicinity of the critical point. The calculations are performed 
with variable /x = for d > 4 and fi = fi^ from Eq. ( [7. 61) (setting e = 1) for d = 3, and 
with variable A taken from Eq. ( |7.5| ) in both cases. The asymptotic behavior for small n is 
presented in the Appendix. The numerical values of B and R in the random case and for 
A^ = 1, presented in Table |, are in full agreement with those obtained in Ref. ||^, while 
those of B for the pure case and A^ = 1 (Table I and Table ||) are in full agreement with 
Ref. ra. 



The random case A^ = 1 for d < 4 can not be considered within the present expansion, 
because of the apparent divergence of /Xq that takes place to the used order in e. Up to now 
there are only numerical values _B = |, (74 = 0.448 and R = §2 = 0.150(7) obtained in p3[ 



through Monte Carlo simulations. What is possible to calculate here are the corresponding 
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values of B and R very close to A^ = 1, e.g. for A^ = 1.001. For completeness these results 
are presented in Table II. More general, one can see that ii N ^ 1, then i? — ;> | and R —>■ 0, 
i.e. the system exhibits SA. This evident discrepancy with the reality is due to the wrong 
assumption that some information about the random case A^ = 1 can be obtained from the 
above formulas in this limiting case. As it was pointed out the correct treatment of the case 
A^ = 1 seems to be a more difficult computational problem. 

The finite size correction to the bulk critical behavior of the cumulants B and R in 
the region Lf^^ = j ^ 1, i.e. away from the critical point, are obtained with the help of 
the asymptotics /i ^ 1, given in the Appendix (Eqs. \K7[ |A^ ). According to the analysis 
presented there, we obtain for Binder's cumulant 



B = 1 



2 

N 



1 + 



3A-1 



3/^2 



+ 



.f^' 



For the cumulant R we get 



R 



/i^ 



+ 



./^' 



^.2) 



The final results can be obtained by replacing A and /i by their respective expressions 
evaluated at the fixed points of the model given in Eqs. (|7.4| ) and ( [7.5|) . So, to the lowest 
order in e, we have 



B= 1 



1 



2 

N 



O 



'e' 



(8.3) 



and 



R 




lA) 



It is interesting to compare Eq. (|8.1| ) with the corresponding result for the pure system 



34 



B. 



pure 




15) 



Up to the lowest order in - they coincide for 1 < A^ < 4, moreover for A^ = 4 — 5 (with 
(5 ^ 1), we have n = jip + 0{t5) and B = -Bpure + 0{e6). The result (^.41) confirms the 



statement that away from the critical point, a strong SA emerges in the system, as i? — > 



with I > 1. 



IX. CONCLUSIONS 

In the present paper we propose a general scheme for the FSS scaling analysis of a finite 
disordered 0{N) system. The method, we use here, is an extension of the field theoretical 
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methods used to analyse FSS properties in pure systems. The nature of the symmetry 
(obtained as a consequence of the use of the rephca trick, which removes the disorder) 
of the model complicates the perturbative structure of the theory in comparison with the 
corresponding 0{N) pure one. Remind that, the final results for the disordered system are 
obtained by making the number of replicas vanishing. Our results concern mainly systems 
with number of components larger than 1 i.e. non Ising systems. Their extension to Ising 
systems requires higher loop calculations, because of the degeneracy of the one loop order 
RG equations. 

Our main results are related to the formulation of the problem for some number of 
components A^ of the fluctuating field for dimensions d > 4 and d = 4 — e. Due to the 
presence of randomness, it is shown that we are dealing with two variables problem with 
scaling variables fi = tL'^/^-u"^/^ and A = A/u. In the mean field regime d > 4 our results 
are a generalization for A^ > 1 of those obtained in for A^ = 1. Evaluating numerically 
the corresponding analytic expressions for the Binder's cumulant B and the relative variance 
R, we demonstrate a monotonic increase of -B as a function of N in both pure and random 
cases and a monotonic decrease of R (to zero for A^ = 4) in the random case (see Table |). 

The e-expansion to first order in e, shows that close to the critical temperature, one can 
express the physical observables, the shift of the critical temperature and the renormalization 
of the coupling constants in terms of n and A. In the random fixed point the parameter 
A takes an universal value (see equation ( [7.5|) ). It is found that the distance, over which 
the bulk critical temperature is shifted, is proportional to L^^^'^^, in agreement with the 
statement of Ref. 0. This result, combined with fenomenological FSS, gives rise to the lack 
of SA. The scaling parameters /i and A, also enter in the final expression for the Binder 
cumulant ( 4.13| ) and the relative variances (^4.14 ), giving explicit expressions for them in the 



different asymptotic regimes. The numerical calculation of the above parameters permits 
also the verification of the FSS and SA (see Table H]). The last is shown to be absent in 
the regime z/r > 2/d, where randomness are relevant and our analysis explicitly shows that 
in this regime the relative variance is always non-zero for the case of stability of the RG 
solutions, i.e. 1 < A^ < 4 and d = 4 — e. 

One can also try to repeat the analysis without the use of the replica-trick. For this aim, 
one needs to define in a proper way the procedure in the zero-mode approximation, when 
rescaling the random part of the Hamiltonian for finite system. One can realize that this is 
the k = part of this term, which will give a contribution to the free-energy. This scheme 
permits formally to end with the same expressions for the shift of the critical temperature 
and the renormalization constants, as we did within the replica formalism. 

In our opinion, the present FSS study can be also applied in the "canonical" case 0, 
where the disorder is characterized by a constant total number of the occupied sites (or 
bonds), instead of the constant average density. We hope that results, similar to the bulk 
case, will also hold in the case of finite geometry, relating in this way our theoretical findings 
with the Monte-Carlo simulations. 
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APPENDIX: FINITE-SIZE SCALING BEHAVIOR OF THE EVEN MOMENTS 

OF THE ORDER PARAMETER 



In this Appendix we present the mathematical details of how to obtain the asymptotic 
of the averages: 



M 



2m 



{M2y 






dx 






x] 



-x-2/2 



where we have introduced the function 



'-N{z) = ON I ""i'^f^ ' exp : 



I^(z) = Sm / (i$$(^-i) 



^$2 + ^$4 
12 



(Ala) 



(Alb) 



(A2) 



The integral in the definition of function Tn{z) given by Eq. (|A2|) , may be evaluated in 
terms of parabolic cylinder functions Dp{z) using the identity [^ 

(A3) 



7 



\S3 



7 



2/?. 



The result is: 



Xn{z) = (1271^)^/4 exp (^ j D^N/2{Vlz). 



Now the above integrand in Eq. ( |Ala| ) can be rewritten in a very simple form 



M2m{x) :-- 



JAr(yU + V^: 



X) 



{n.'y 



D. 



■m-N/2 



((/i + v^a;)\/3) 



D_M/2{{^ + y/Xx)V^) 



(A4) 



(A5) 



For small // <^ 1 (i.e. in the vicinity of the critical point) the asymptotic form of the 
ratio ( |A5|) is given by 



M2m{x) = (l27r2) 



f jL'_„_Ar/2(xV3A) 



D_Ar/2(a;v3A) 



iV- 



+0{iif}. 



D-N/2-lix\/3X)D_rn-N/2ixV3X) 



(2m + A^) 



D_m-Ar/2~i(a;v3A) 



^-7v/2(a;v3A) 



(A6) 
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In the mean-field regime and at the critical point we have /i = and Ai2m is equal to 
the first term in the r.h.s. of Eq. (|A6|) . 

For large /i S> 1, the asymptotic behavior of the ratio Ai2m is obtained with the help of 
the well known Watson's Lemma: 



Lemma: (see for example [Q) 

Suppose a > 0, f3 > and f{x) an analytic function in a neighborhood of x = 0, 

f{x) = £ ^^x\ \x\ < R. 



k=0 



k\ 



and that 



\f{x)\<c,e' 
for positive constants Ci, c^,. Then: 



X 



X' 



^-ie-^"7(x)dx ~ - E ^ 



-| oo 

a 



k=0 



xe [R,X], 

[a k\ ' 



as s -^ oo in the sector \ args| < f . 



According to this Lemma, from Eqs. ( |A^ ) (with z = fj. + xvA) and ( |A5|) , we have 



M2(x) = (127r2 



1/2 N 



and 



MAx) = 12tt' 



,N{N + 2) 



/^^ 



_ xV\ Gx^A -N-2 /I 

2x^A Ox^A - A^ - 3 



(A7) 



1- 



fJ' 



+ 



3/i2 



O 



./"' 



(A8) 



Using the asymptotic of A^2 and M^ for large /z, we can get the behavior of the cumulants 
R and B in both cases d > 4 and d = 4 — e. They are given by 



3 V A^ 



1 



3A-1 

3/i2 



<i.: 



(A9) 



For the cumulant R we get 



/i2 l/i3 



(AlO) 
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TABLES 



TABLE L Numerical values for the Binder cumulant B from Eq. ( [4.13 ) and the relative 
variance R from Eq. ([4.14|) in the mean-field regime i.e. d > 4. 



Random 



Pure 



N 



B 



R 



B 



R 



1 

2 
3 

4 



0.216368 
0.451486 
0.533513 
0.575587 



0.310240 

0.111381 

0.038365 





0.270520 
0.476401 
0.543053 
0.575587 










TABLE II. Numerical values for the Binder cumulant B from Eq. (4.13) and the relative 
variance R from Eq. ( [4.14| ) at d = 3 



Random 



Pure 



N 



B 



R 



B 



R 



1 
1.001 
2 
3 
4 



0.666334 
0.602793 
0.625783 
0.640628 



0.000427 

0.061279 

0.022688 





0.400024 
0.400328 
0.547496 
0.592813 
0.614002 
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FIGURES 




FIG. 1. One loop contributions to the reduced temperature t 
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FIG. 2. One loop contributions to the couplings u and A 
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